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abstract 


A  formula  is  derived  for  the  force  between  two  parallel  coaxial  coils  whose 
winding  channels  have  rectangular  cross  sections.  The  expression  for  the  maxi- 
mum force  also  is  obtained.  The  linear  dimensions  of  the  cross  section  of  the 
windings  are  all  small  compared  with  the  mean  radii  of  the  coils  and  sixth-order 
terms  in  these  ratios  are  neglected.  The  principal  term  in  the  force  corresponds 
to  two  linear  circular  currents  in  the  mean  positions  of  the  coils.  The  remaining 
terms  take  account  of  (a)  the  finite  cross  section  of  the  winding  channels,  (6)  the 
discrete  nature  of  the  windings  (round  insulated  wires)  with  waste  space  between 
them,  and  (c)  the  nonuniformity  of  current  over  the  section  of  the  wire.  The 
effects  (&)  and  (c)  are  eliminated  when  the  ratio  of  the  mean  radii  of  the  coils  is 
determined  electrically. 
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I.  INTRODUCTION 

The  absolute  measurement  of  an  electric  current  by  the  use  of  a 
certain  type  of  current  balance  requires  the  computation  of  the  force  g 
between  two  coaxial  coils  in  which  electric  currents  Ix  and  I2  are  circu- 
lating, and  when  the  coils  are  so  spaced  as  to  make  this  force  a  maxi- 
mum, which  is  then  indicated  by  %m. 

If  the  distance  between  their  mean  planes  is  in  the  general  case  z 
(and  zm  when  the  force  is  a  maximum),  and  if  their  mean  radii  are 
di  and  a2,  and  if  the  axial  breadth  26 *  and  radial  depth  2ct  (for  ^  =  1,2) 
of  the  rectangular  cross  sections  of  the  winding  channel  are  such  that 
the  quantities  bt  and  ct  are  small  compared  to  either  of  the  three  finite 
quantities  a1}  a2,  or  z,  then  the  principal  part  of  this  force  is  ni,  n2,  I\ 
I2  F  where  nt  is  the  number  of  turns  and  F  is  the  force  between  two 
circular,  coaxial,  unit-current  filaments  having  the  mean  radii  and 
lying  in  the  mean  planes  Z\  and  z2  of  the  coils.  This  force  i^is  a  func- 
tion of  two  dimensionless  variables  x  and  a. 

a= J  <1,  '*-£,  A*=aZ+aS  (1) 
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where  z  =  zx  —  z2.  When  x  and  a  are  known,  F  (x,  a)  may  be  computed 
with  great  precision  by  Maxwell's  formula  in  terms  of  the  complete 
elliptic  integrals  K  and  E  with  modulus  k  where 


„           4«,«2                2V1-02        ^nrtfl_l-«2 

(2) 

'"      {<h  +  <hY+z>     l+a2+Vl-^         "      l+«2 

F(x,a)  =  ^Jj     k¥E    2K]-2V2A(l-j32)-l^'(F) 

(3) 

where  4>'(P)  denotes  d4>ldl<?  and 

♦-a^-^~ 

(4) 

The  function  <£  satisfies  the  differential  equation 

V{l-W)4>'' =¥<}>'  +  \<t>  (5) 

The  force  3  between  the  actual  coils  obtained  by  integration  over 
their  cross  sections   will  be  expressible  in  the  form 

3(*,a)  =n1n2iiI2F(x,a)[l  +  A2(x,a)  +  A4(z,a<)  .  .  .]  (6) 

where  the  functions  of  x  and  a  designated  by  A2  and  A4  are  correction 
terms  of  the  second  and  fourth  order  in  the  small  constant  ratios 
bi/ai,  cjdi,  etc.,  and  are  to  be  found. 

II.  THE  FORCE  WHEN  THE  CURRENT  UNIFORMLY  FILLS 
THE  WINDING  SPACES 

The  force  between  two  current  filaments  of  radii  rx  and  r2  planes 
z[  and  z2  may  be  expanded  by  Taylor's  theorem  in  terms  of  the 
small  quantities  z'i~z1}  z'2  —  z2,  ti^-di,  r2  —  a2.  Integrating  this  over 
the  two  cross  sections  gives  on  neglecting  terms  of  higher  order  than 
the  fourth, 


n 


^jj2  =  F+  J  j(6J  +  bl)D>3F+  cWa\F+  cJZ%*} 

+  m{l3{ht  +  bt)  +  ^2M]DtF+S(ciD^F+ciDa\F)    (7) 

+  10clc£DIlDl2F+  10(bl  +  bl)(clD2aiDlF+cm2DlF)\ 

where  the  partial  derivatives  of  F  imply  it  is  a  function  of  three 
variables  z,  au  and  a2. 

A  first  step  in  reducing  this  to  a  workable  formula  is  to  express 
all  the  partial  derivatives  in  terms  of  ^-derivatives  and  then  these  in 
terms  of  ^-derivatives.     Let  Fl,  F11,  etc.,  denote  DXF,  DX2F,  etc.,  so 
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that  DznF=Fw/An.     Then  it  is  found  by  the  use  of  equations  (1)  to 
(5)  that 


(8) 


By  use  of  these,  all  the  terms  of  second  order  in  (7)  may  be  expressed 
in  terms  of  F,  F1,  and  F11. 

For  the  partial  derivatives  of  fourth  order  which  appear  in  (7),  one 
finds. 


712  7}  27^=  _± 

Ua2J-Jz  £  —  a  i 


(9) 


A'      2A 
It  is  also  found  that 


mxF= 


~AT  +  'M?l 


2DitD2tF 


and 


Wiv       jpii 


(10) 


Bixm2F= 


IF 


DilD>F-DlpiE+2^l(l+i*)F11  +  ZxF*\    (11) 


By  use  of  (9)  to  (11)  the  term  of  fourth  order  in  (7)  is  expressible  in 
terms  of  F,  F1,  F11  .  .  .  Flv.  It  is  then  possible  to  express  F11, 
F111,  and  FlY  in  terms  of  F  and  F1,  for  it  is  found  from  (1)  to  (5) 
that  if  ax  and  a2  are  constant  the  force  F  satisfies  the  ordinary  differ- 
ential equation  in  x. 

(xi  +  2x2+p2)x2FII+(5xi  +  2x2-Sp2)xFI-(2x2-S(32)F=0      (11)' 

From  this  and  the  equation  found  by  differentiating  it  successively, 
all  derivatives  of  higher  order  than  the  first  may  be  expressed  in 
terms  of  F  and  F1.  However,  for  the  purpose  of  computation,  it  is 
simpler  not  to  make  this  elimination  but  to  compute  in  succession 
the  functions  of  x  and  a,  \u  X2  .  .  .  X4,  which  are  defined  by 


Xi^ 


X2  = 


xF1 


x2F11  J6$2-2x2+  (5a;4  +  2a;2  -  3  ff2)  \x 
F  x*  +  2x2+p2 

[ii 


_xzFlu  _  4a;2  -  (25a;4  +  4a;2)  \  +  (1  la;4  +  10a;2  -  ff2)  X2 
A3~     F     "  a;4  +  2a:2+/32 


X4^ 


F 

x*FIV 


a;2[75a;2X1-3(23a;2  +  8)X2+16(a;2+l)X3] 
a;4  +  2a;2+j82 


(12) 
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In  this  way  the  expression  (7)  leads  to  the  form  (6)  where 

*<ft  «)=li2j-[5-|>-[(^-«|+(x2+»|>1 


2{b\-c\  +  b22-c 


lA*  ZJX2)     (13) 

A4(z,  «)  -ggJ^JB^-  (S2x2-j3B1)X3-[5(^)V+53^-3^I]x2 

+  3[<^)*  +  20»>>-W»>l     (14) 
where 

5i  - 1  [5  (6?  +  b\  -  el)  -  3c?]  - 1  [5  (6?  +  b\  -  c?)  -  3d] 

B2=|[5(6?  +  61-cI)-3c?]  +  |[56?+6I-Cf)-3c!]  (15) 

#3=4!  [9cf+15cI-10(6?  +  6I)]+|  [9cS  +  15cf-10(6?  +  6|)] 
a\  a% 

B4-2;2[-3(6}  +  c}  +  6J  +  cJ)  +  10(6M  +  6^)-10(6;-cD(6I-cI)] 
III.  THE  MAXIMUM  FORCE 

In  the  absolute  measurement  of  current,  the  distance  between  the 
mean  planes  of  the  coils  is  adjusted  by  trial  to  a  value  zm  which  makes 
%  a  maximum.  The  corresponding  value  of  x  is  xm.  The  value  z0 
(and  x0  corresponding)  which  makes  g  a  maximum  is  slightly  different 
on  account  of  the  correction  terms  A2  and  A4.  If  primes  denote 
^-derivatives  as  before,  the  equation  to  determine  xm  is 

ffjxm,  a)  =81(^>  g)  i  ^(Xm,  a)  +  A41(a;m,  a:)^0  ," 

g(xm,  a)       g(zOT,  a)  1  +  A2(zm,  a) 

whereas  the  corresponding  equation  determining  x0  is 


F(z0a) 


=  0  (16)' 


Writing  Xm=xo  +  8x0  in  (16)  and  expanding  the  functions  on  the  right 
in  ascending  powers  of  8x0}  and  taking  account  of  (16)',  shows  that 
8x0  is  a  second-order  quantity  whose  principal  part  is 

8Xo^-^M  (17) 

Consequently  the  maximum  force  gTO  =  g  (xm)a)  is  given  to  the  fourth 
order  inclusive  by 

^m  =  n1nJ1I2F(x0,a)^l  +  A2(x0)a)  +  Ai(x0)a)  +  ^ffiffi*]     (18) 
where 
x0A2'(x0,a)  =I^-2|2^(|-|)  +  [«-«  |+  (4+/3)  |]  X2(x0)«) 
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+*<*-*+*-*  ^  (19) 

*<»,«)-  T^j-(g-|>-2  W-^+«-^  ^      (20) 

The  term  A4(£0,a)  is  computed  from  (14)  and  (15)  for  the  value  x  =  x0 
which  makes  Xi  =  0  so  that  by  (12)  the  functions  X2,  X3,  and  X4  become 
rational  fractions  in  x02  and  |82. 

IV.  THE  FRACTIONAL  ERROR  IN  THE  MAXIMUM 

FORCE 

In  evaluating  the  small  terms  A2,  etc.,  the  axial  width  2bt  and  radial 
depth  2ct  of  the  winding  section  must  be  measured.  The  mean  radii 
of  the  coils  (&i  and  a2)  may  also  be  obtained  from  measurements  with 
sufficient  precision  for  this  part  of  the  computation.  For  the  prin- 
cipal term  F  (x0)  these  enter  only  through  their  ratio  a  =  a2/ar  and  this 
must  be  known  with  the  highest  precision  although  x0  need  be  known 
only  with  moderate  exactness  since  its  error  makes  an  uncertainty  of 
second  order  in  the  force.  The  fractional  error  in  %m  due  to  this 
source  is 

(f ).-4  ?Sfe©'— <£?«--■'  w 

approximately  so  that  an  error  of  0.1  percent  in  x0  would  only  make  an 
error  of  1.5  parts  in  a  million  in  the  force  gw.  It  is  evident,  therefore, 
that  when  a  is  known,  x0  may  be  obtained  with  ample  precision  by 
interpolation  in  the  table  4  given  by  Dr.  Grover.1  This  table  gives 
ym=Zo/a!SO  that  by  (1) 

*-=v&  (22) 

An  estimate  of  the  fractional  error  introduced  into  the  force  (equation 
(18))  by  all  errors  in  measurement  of  the  dimensions  of  the  coils  is 
found  by  taking  the  first-order  variation  of  this  equation.  In  doing 
this  the  8x0  is  omitted  since  by  (21)  it  produces  no  variation  of  first 
order.  Also  the  terms  of  fourth  order  in  (18)  may  obviously  be 
unvaried  (since  they  are  almost  negligible  themselves).  Further- 
more, in  the  term  A2  the  errors  8b  t  and  8c <  in  the  measurement  of  the 
cross-sectional  dimensions  of  the  coils  will  be  dominant,  while  those 
caused  by  uncertainty  in  ai}  a2,  and  x0  will  be  here  of  no  practical 
importance. 

With  this  understanding,  the  first-order  variation  in  %m  is 


(1  +  j8)  X2  b\  hbx     (1  -  0) X2  hi  8h] 

3        a\  &i  3        a\  b2\ 


(23) 


i  F.  W.  Grover:  Bui.  Bu.  Stand,  vol.  12,  Sci.  Paper  No.  255,  p.  372;  1915. 
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This  is  not  final ;  it  is  necessary  to  express  8a/ a  in  terms  of  quantities 
more  directly  measurable.  The  ratio  a  is  obtained  by  measuring  the 
current  ratio  T2\TX  when  the  coils  carrying  these  currents  are  made 
coplanar  and  produce  equal  axial  components  of  magnetic  field  H  at 
their  common  center. 

P  If  H  be  computed  on  the  same  assumption  which  has  been  made  in 
deriving  the  formula  for  $m,  which  is  that  the  current  uniformly  fills 
the  winding  channels,  it  is  found  that 

with  a  similar  expression  for  H2,  where 

J*»f -f  and  ^  =  ^-6?cf  +  f6?(f=l,2)  (24)a 

When  Hx  =  H2  the  ratio  a  is  given  (to  the  fourth  order)  by 

-*{'-'J(*-^D«-«^a?*]|    » 

where  r  is  the  measured  ratio  n2/2'Mi^i-  The  first  variation  of  a  is 
found  from  this  (with  the  same  assumptions  as  were  used  in  getting 
(23)  ). 

8a  _  8r     8Ji     8J2  _ 

a      t      a\       a\ 

r     3  \    a\  d         a\  c2      a\  hx      a\  b2\  ^    ) 

Introducing  this  value  of  8a/ a  into  (23)  gives  finally  the  fractional 
error  in  gm 


8 


[cj8^_M8b2l] 
[aic2     a22b2J\ 


(27) 


This  expression  is  symmetric  with  respect  to  the  two  coils.  This  may 
be  seen  by  noting  that  $  and  6tJt  are  antisymmetric,  while  xl  and 
X2  are  symmetric  with  respect  to  the  subscripts  1  and  2. 

V.  CORRECTION  FOR  CURRENT  DISTRIBUTION 

In  this  section  are  found  corrections  to  gTO  and  also  to  a  (in  both 
cases  of  fourth  order),  which  are  made  necessary  by  the  facts  that  (a) 
the  wires  are  round,  leaving  some  waste  space  in  the  winding  chan- 
nels, and  (b)  the  current  density  may  not  be  uniformly  distributed 
over  the  section  of  each  wire,  as  is  generally  assumed  for  steady  cur- 
rents. It  is  in  fact  possible  that  a  theoretical  examination  might 
lend  more  probability  to  the  so-called  natural  distribution  in  which 
the  current  density  at  any  point  in  the  wire  is  inversely  proportional 
to  its  distance  from  the  axis  of  the  turn. 
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It  will  presently  be  seen  that  the  correction  due  to  current  distri- 
bution which  must  be  made  in  %m  is  compensated  for  by  the  correction 
which  must  be  made  in  a  when  it  is  determined  electrically — the 
resultant  of  the  two  corrections  is  of  higher  order  than  the  fourth. 

To  obtain  the  correction  (a)  imagine  each  coil  to  consist  of  nt  turns 
of  square  wire  all  in  contact,  each  having  a  uniform  current  distri- 
bution. Let  gs  denote  the  force  between  two  such  turns,  one  in  coil 
no.  1,  the  other  in  no.  2.  Since  there  is  no  waste  space,  the  force  g, 
equation  (6),  is  2g„  the  summation  being  extended  to  all  possible 
pairs  of  turns,  one  of  which  is  in  no.  1,  the  other  in  no.  2. 

The  force  S  +  Sg1  between  the  actual  coils  which  consist  of  discrete 
windings  with  wire  of  radius  pt  will  be  the  corresponding  sum 

where  gc  is  the  force  between  the  two  turns  of  circular  wire  which 
have  the  same  mean  plane  and  same  mean  radius  as  the  pair  of  turns 
of  square  wire  to  which  gs  refers,  the  current  being  assumed  to  be 
uniform  in  both  cases.  If  their  mean  planes  are  Z\  and  z2  and  mean 
radii  are  rx  and  r2,  it  is  found  by  an  expansion  similar  to  (7)  that 

^g^(($'-|^>+x{(x-«^/+(^+«3^]j    (28) 

+    terms  of  higher  order   than  the  fourth  (where   x=  1  A  2  =-j\ 

The  term  retained  is  in  fact  of  the  fourth  order  since  the  small  quan- 
ties  e/  are  of  the  second  order,  defined  by 


4}-%&)  (28,) 


a$mbi 

r\   nc 

where  2+  /^£l  is  the  side  of  the  square  wire  section  and  2  pt  the  dia- 

meter  of  the  circular  wire.     Since  irplni  is  the  total   metallic  cross 

section,  the  fraction  1  —  jr   i  is  approximately  the  fraction  of  cross 

section  which  is  waste  space.     For  ordinarily  close  windings,  this 

will  be  of  the  order  of  magnitude  (but  never  less  than)  1  — j  =  0.21. 

Hence  if  the  ratio  of  the  diameter  of  the  wire  to  a  side  of  the  cross 
section  of  the  winding  channel  is  of  the  same  order  as  bi/ai,  i.e.,  of  the 
first  order,  then  1/%  and  e'  will  be  of  the  second  order  and  the  terms 
above  representing  gc—  gs  are  of  the  fourth  order,  and  we  may 
place  ai/r2i  =  l,  x  =  x0,  which  gives 


«-!&-£;  (i=1'2)  (29) 


and 


%<=%-IihFi2xi a)  KaH'^h) (since  h(x°' a)=o)' 
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Summing  the  nxn2  elementary  forces  between  pairs  of  turns  of  this 
type  gives 

Ztic  —  tfm^ofim  -Zfis  ^2 ^\a2  1~a2  2) 

Hence,  since  2gs  =  gTO,  we  find  for  this  fourth-order  correction 


OjSrn 

5w 


iii(«?el~3f0  (30) 


This  takes  account  of  the  fact  that  the  coils  are  wound  with  wire  of 
circular  section  having  a  certain  amount  of  waste  space  between 
them.  This  8^'m  must  be  added  to  the  second  member  of  (18)  to  get 
the  maximum  force  between  the  actual  coils. 

For  the  correction  (b),  consider  the  current  density  %i  in  a  round 
wire  in  coil  no.  1  to  be  a  function  of  the  distance  r/  of  a  point  in  this 
turn  from  the  axis  of  the  turn.  Then  expanding  to  the  second  order 
we  find 

For  uniform  distribution  (7i  =  (7'i  =  0;  for  the  natural  distribution 
C\  =  Ci  =  —  1.  If  Fs'  is  the  force  between  two  circular  turns  when  the 
current  distribution  is  given  by  (31)  and  Fs  as  defined  above,  the  force 
between  them  when  their  current  is  uniformly  distributed,  we  find 
with  an  expansion  similar  to  (7)  by  integrating  over  their  sections 
that  Fs'  —  Fs  is  a  small  quantity  of  the  fourth  order  which  is  obtained 
by  replacing  the  e/  in  the  second  member  of  (28)  by  e-  where 

e';  =  ftf,|  (32) 

Hence  by  a  similar  argument  the  correction  to  the  force  Fm  due  to 
this  cause  is  given  by 


sf:        (3 


Fm  12x1 


6h-oH)  <33) 


The  correction  5*gTO  =  5gJ,  +  5g^  due  to  waste  space  and  to  non- 
uniformity  of  current  in  each  wire  is  given  by 

^=-T^Xrtei~^v  where  ei=ei+€i-        (34) 

The  details  in  computing  these  two  effects  on  the^magnetic  field  H 
produced  by  the  coils  are  so  similar  to  those  just  sketched  that  they 
need  not  be  given.  The  final  result  is  that  (24)  is  only  altered  by  a 
change  in  jt  so  that  (24  )a  has  to  be  replaced  by 


■'.-K'-sH 


(35) 
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The  result  of  this  change  is  to  add  a  term  5* a/ a  to  the  second  member 
of  (26)  where 

5*a     \(c\        c\    \ 

The  contribution  of  S* a/ a  to  the  second  member  of  (23)  is 


12z; 


m«-m 


which  is  exactly  canceled  by  the  correction  (34)  to  gTO  which  must 
be  made  at  the  same  time.  Thus,  after  taking  account  of  these 
two  distribution  effects  (to  the  fourth  order),  we  come  back  to  pre- 
cisely the  earlier  equation  (27)  for  the  fractional  error  in  the 
force.  But  this  equation  (27)  was  derived  from  the  expression  (18) 
for  the  maximum  force.  Moreover  (18)  was  derived  by  assuming 
that  the  current  uniformly  filled  the  winding  channels,  that  is,  by 
neglecting  both  of  these  distribution  effects.  We  arrive  then  at  the 
following  conclusion. 

VI.  CONCLUSION 

If  the  ratio  a  of  the  mean  radii  of  the  coils  of  the  current  balance 
is  determined  electrically  and  this  value  is  used  in  the  computation 
of  the  maximum  force  between  the  coils  using  equation  (18),  the 
force  thus  computed  will  differ  from  the  maximum  force  between 
the  actual  coils  by  an  amount  which,  compared  to  this  force,  is  a 
small  quantity  of  order  higher  than  the  fourth  in  the  small  quantities 

— »  —  i  --»  and  — »  where  2  bt  and  2  c{  are  the  axial  breadth  and  radial 
«i  a\  a2  a2 

depth  respectively  of  the  rectangular  cross  sections  of  the  winding 

channels — provided  also  that  the  number  of  turns  n*  in  each  coil  is 

so  large  that  1/nt  is  in  each  case  i  =  l  and  t  =  2  a  small  quantity  of 

the  order  of  (bi/ai)2  or  smaller.     These  conditions  are  met  in  practice. 

The  effects  of  errors  in  measurement  of  cross  sectional  dimensions 

may  be  computed  by  equation  (27).     This  conclusion  rests  upon  the 

derivation  of  equation  (18)  and  also  upon  the  investigation  of  section 

V,  where  account  is  taken  of  the  facts  that  the  actual  coils  consist 

of  many  turns  of  round  wire  with  waste  space  between  them,  and 

that  the  current  may  not  be  uniformly  distributed  over  the  section 

of  the  wires  (although  the  distribution  is  assumed  to  be  of  the  same 

nature  in   all  wires).     This  section   appears,   therefore,   to   give   a 

greater  validity  to  the  force  equation  (18)  than  was  to  be  expected 

from  the  mode  of  its  derivation  for  it  was  derived  on  the  assumption 

that  the  current  uniformly  filled  the  winding  channels.     The  crux 

of  the  matter  is  the  electrical  determination  of  the  ratio  of  the  radii. 

If  this  ratio   a  were  determinable  with   the   desired  precision  by 

geometrical  measurements  of  the  coils,  the  conclusion  stated  above 

would  not  be  true.     In  that  case,  the  computation  of  the  force  would 

require  a  knowledge  of  the  correction   term   5*Fm  of   (34),  which 

would  then  be  added  to  the  second  member  of  (18). 

Washington,  August  8,  1933. 


